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■^j- | Knizhnik-Zamolodchikov-Bernard equations for twisted conformal blocks on compact Riemann sur- 

■ faces with marked points are written explicitly in a general projective structure in terms of corre- 

0^ | lation functions in the theory of twisted b-c systems. It is checked that on the moduli space the 

equations provide a flat connection with the spectral parameter. 

O 

Q 

0^ ' In addition to the conformal symmetry, the Wess-Zumino-Novikov-Witten (WZNW) model possesses 

(N ■ the symmetry of an afflne Lie algebra Q [1]. The Virasoro algebra is embedded in U(Q) by the Sugawara 

' construction. This additional symmetry leads to certain equations for the conformal blocks in the WZNW 

theory. When the theory is defined on the sphere, these are the well-known Knizhnik-Zamolodchikov (KZ) 
equations [2]: 

§i ijk + *TF S ^7=^) <*i(*) • • • <M*0> = o. <i.i, 



c3 



1. Introduction. 



These equations relate the dependence of the conformal block ($ 1 (z 1 ) . . . $„(z„)) on the positions of the 

(1.2) 



| marked points z a to the action of the currents j a (z) on the fields <& a ( z a 

Jw * m 

We may also speak of the KZ equations as of the connection 



V - 9 I 1 V af> (13) 

u ; 

on the bundle of conformal blocks over the moduli space of the sphere with marked points. 

The subject of this paper is the generalization of the KZ equations (1.1) to surfaces of higher genera 
and the discussion of their properties. First these equations for the surfaces of nonzero genera were obtained 
by Bernard [3,4]. To define the action of the zero modes of the current he included twists on the handles 
of the surface. Conformal blocks become functions of the twists, and their derivatives along the twists give 
the action of the zero modes of the current. The moduli space contains not only the positions of the marked 
points, but also the moduli of the complex structure of the surface itself. The connection form becomes a 
differential operator with respect to the twists. We shall call such a system of the equations the Knizhnik- 
Zamolodchikov-Bernard (KZB) equations, and the corresponding connection — the KZB connection. 

In the paper we write out the KZB equations in a simple form. All terms in the equations can be 
expressed explicitly as Poincare series in the Schottky parametrization, and at the same time admit an in- 
variant description in an arbitrary projective structure. We also briefly discuss the transformation properties 
of the equations as the projective structure changes. A remarkable feature of the KZB equations is their 
relationship to the twisted b-c system of spin 1. After multiplying the conformal block by the square root 
of the b-c holomorphic partition function, the connection depends on the level k only through a spectral 
parameter (fc + ft*) -1 ; besides the connection form becomes symmetric. It was pointed out by Losev that 
this relationship to the b-c systems can be explained in terms of the BRST construction in the G/G coset 
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model [5] . The presence of the spectral parameter in the connection ensures that from its flatness at integer 
k (for "physical" reasons) the flatness at an arbitrary k will follow [5] . This imposes certain strong conditions 
on the connection form. We explicitly check these conditions and prove the flatness. 

At higher genera the KZB equations contain an interesting "potential" term. This term vanishes at 
genus one (on a torus), and Bernard originally claimed that it was zero at any genus [4]. However there exist 
indications that the potential should not vanish at higher genera [5]. So far we cannot say much about the 
potential term, except its closeness as a 1-form on the moduli space, which is one of the conditions for the 
compatibility of the KZB equations (cf. [8]). 

The paper is organized as follows. In section 2 we introduce the twisted WZNW model on a compact 
Riemann surface. In section 3 we recall how the stress-energy tensor defines the connection on the moduli 
space and how this connection depends on the choice of the projective structure. Section 4 deals with twisted 
1-forms on the surface; we also define the twisted b-c system of spin 1 in this section. These sections provide 
a necessary kit for working with the KZB equations. Here we tried to follow the ideas and the notation of 
Bernard [3,4] and Losev [5]. Finally, section 5 contains the main results of the paper. In this section we write 
out the equations and discuss their properties. Section 6 summarizes the discussion and presents several 
questions for further investigation. Some auxiliary information and calculations are gathered in Appendices. 
The Schottky parametrization of Riemann surfaces is reviewed in Appendix A. In Appendix B we derive the 
KZB equations. Appendix C contains the proof that the potential term is a closed 1-form on the moduli 
space. In Appendix D we commute the derivatives with respect to the positions of marked points. In 
Appendix E we check the compatibility of the equations for the partition function which yields the proof of 
the compatibility in the general case. 

2. Twisted WZNW theory. 

We consider the WZNW model at level k on a compact Riemann surface of genus N . The model is 
defined for a simple compact Lie group G with the Lie algebra Q [2] . We shall deal only with the holomorphic 
part of the theory, thus all fields are holomorphic except at locations of other fields. The holomorphic theory 
contains the currents j(z) taking values in the Lie algebra Q, the stress-energy tensor T(z), and primary 
fields &(z). The primary fields take values in irreducible finite dimensional representations of Q and are 
multivalued on the surface. The fields obey the operator product expansion (OPE): 

f(z)j\w) = -kj-—^ j^fM + 0(1), (2.1) 

f(z)<i>(w) = j^t a $(w) + 0(1), (2.2) 

where S ab is the invariant bilinear form on Q, t a is the action of a £ Q in the representation of the primary 
field <&, f abc are the structure constants of Q. We shall work in an orthonormal basis of Q normalized so that 

jabejabd = 2h*8 cd , (2.3) 

where h* is the dual Coxeter number of Q. The stress-energy tensor T(z) is expressed in terms of the currents 
by the Sugawara construction: 

T(z) = „ 1 , . : j a (z)j a (z) : = 1 - - lim ( j a (z)j a (w) + kdlmG '' \ . ( 2 .4) 

v ' 2(k + h*) J y IJ y 1 2(k + h*)w^z\ J y /J y ' (z - w) 2 J y ' 

To be specific about notation, we shall label normalized correlation functions (divided by the partition 
function) with the subscript N, like (X) N , non-normalized correlators being without any subscript: (X) = 
(X)mZ, where Z = (1) is the partition function. 

Our ultimate goal is to study the connection on the bundle of correlation functions over the moduli 
space. To write this connection is the same as to compute the correlation functions with the insertion of 
T(z) (see section 3). Since T(z) is constructed of the currents, to write the KZB connection we first need to 
compute correlation functions with the currents inserted. On the sphere any correlation function (j a (z)X) 
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is a 1-form in z which can be reconstructed from its singularities (singular terms in Laurent expansion at 
the poles). On Riemann surfaces of higher genera there exist global holomorphic 1-forms, and to restore 
a mcromorphic 1-form, we need extra information, besides singularities. For example, any meromorphic 
1-form is uniquely defined by its singularities and integrals over A-cycles. Here we encounter difficulties, 
since zero modes 

£(j a (z)X)dz (2.5) 

are not determined by the OPE. For this reason Bernard suggested to include twists on A-cycles [3,4]. Given 
a set of ./V elements of G (N is the genus of the surface) 

.g^expQT^), i=l,...,N, (2.6) 

a 

we insert in all correlation functions the twists 

gi=expf &j a {z)dz. (2.7) 

From now on we shall suppress the twists in notation, thus the correlators (X n£Li 9i) wm be written simply 
as {X). 

Following Bernard, we write equations not for a single partition function or a correlation function, but 
for them as functions of the twists. This allows us to express zero modes of the currents as the derivatives 
along the twists: 

U dzf{z)X)^C a {X) 1 (2.8) 
where C la is the right-invariant derivative along the i-th twist: 

J(g ll ...,^ ta g ll ...,g N ) (2.9) 

for any function / of N group elements gi. 

3. 1-forms on the moduli space and projective structures 
on the surface. 



jC ia f( 9l ,...,g N ) = — 



The variation of a correlation function under an infinitesimal shift of the moduli is described by inserting 
the stress-energy tensor T(z). Indeed, let us consider an infinitesimal change of complex structure induced 
by a coordinate transformation 

z i— » z + e(z, z), (3.1) 
where e(z,z) is defined only locally (otherwise it would give a reparametrization) . The Beltrami differential 

[i(z, z) = Bs(z, z) (3.2) 

is defined globally and has the transformation properties of a (-l,l)-form on the surface. At the same time it 
may be thought of as a tangent vector to the moduli space. To compute the variation of a (non-normalized) 
correlation function under the change of moduli described by /i, we need to insert the stress-energy tensor 
coupled to fj, inside the correlator: 

S„(X}= [ (T(z)X) fi(z,z)dzdz. (3.3) 

This expression would be well-defined if T(z) were a 2-differential on the surface E. In fact, the holomorphic 
stress-energy tensor T(z) is not a 2-differential, but transforms with the Schwarzian term: 

TH=(^) 2 T(z) + ^{z; w} , (3.4) 
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where 

, fd 3 z /dz\ 3 / d 2 z idz\ 2 
{z;w}= bLj-2W^J ' (3 ' 5) 

c is the Virasoro central charge. The Schwarzian derivative {z; w} vanishes for projective transformations, 
therefore after fixing the projective structure T(z) becomes a 2-differential. Thus the values of the correlation 
function depend on the choice of the family of projective structures on the surfaces [7]. Physically, this is due 
to the anomaly; prescriptionally, this dependence appears in the regularization of the stress-energy tensors 
(2.4), (4.18). One easily checks that the normal ordering in (2.4) and (4.18) depends on the choice of the 
local coordinate, and that this dependence reproduces exactly the transformation law (3.4). 

The projective structures should be chosen in such a way that the system of equations (3.3) is compatible. 
Then the KZB equations would be compatible, i.e. they would define a flat connection (not just projectively 
flat). Such families of projective structures exist. One of them is the projective structure defined by the 
Schottky paramctrization (see Appendix A). The Schottky parametrization is also convenient for writing 
out explicit formulas, and we shall present the expressions for all terms of the equations in the Schottky 
representation. However, we must stress that the whole treatment is paramctrization independent, and can 
be performed in any "compatible" projective structure (i.e. such that the equations (3.3) are integrable). 
Any change of the projective structure is defined by adding a 1-form on the moduli space to T(z), therefore 
one can easily construct a "non-compatible" projective structure by adding a non-closed 1-form. From now 
on we fix a "compatible" projective structure and speak of T(z) as of a 2-differential, and of correlation 
functions as of functions on the moduli space. 

Remark now that T(z) being a 2-differential on the surface serves as a 1-form on the moduli space [7]. 
Indeed, by (3.3) it defines the response to any change of the moduli. In particular, if we change the moduli 
by moving a marked point £, then 

f<- Y > = |li< T < 2) - Y >- <"> 

It allows us to view any meromorphic 2-diffcrcntial with poles at marked points as a 1-form on the moduli 
space. For a surface of genus N > 2 with n marked points there exist (3 AT — 3 + n) linearly independent 
meromorphic 2-differentials with simple poles at marked points. They form a basis in the cotangent space 
to the moduli space of the surface with marked points. Remark that the second-order residues of the stress- 
energy tensor are fixed by the field dimensions, higher orders vanish if the fields are primary with respect to 
Virasoro algebra. 

Due to the equivalence between 2-differcntials on the surface and 1-forms on the moduli space, the 
differential d m on the moduli space mapping functions to 1-forms can be treated as an operator mapping 
functions on the moduli space to 2-differentials on the surface. We denote this operator by d m (z) so that 

S^F — J n(z, z) d m {z)F dzdz (3-7) 

for any function F of the moduli. Notice that for F being a correlation function, 

d m (z)(X) = (T(z)X). (3.8) 



4. Twisted 1-forms on the surface and 
the twisted b-c system. 

Twisted 1-forms. 

Any correlation function of the form (j a (z)X) is a meromorphic twisted 1-form in z with values in the 
Lie algebra Q (we shall also call it a twisted 1-form in z,a). Twisting means that passing across A-cycles 
(where the twists gn are placed) this 1-form is conjugated by the corresponding twists. More formally twisted 
1-forms can be described with the help of the "Schottky-type" covering. Namely, let us consider the covering 
S* of the surface £ such that the set of B-generators of 7Ti(£) lifts to an infinite tree on £*, and A-generators 
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lift again to closed loops. Recall that the fundamental group 7Ti(E) is generated by the 2N elements Ai and 
Bi, i = 1, . . . , N, with the only defining relation 

A 1 B 1 A^ 1 B^A 2 B 2 A^ 1 B^ 1 . . . A N B N A^B^ = 1. (4.1) 

7i"i(E*) is the subgroup of £ generated by all Ai and group elements conjugate to Ai. Evidently, it is a 
normal subgroup, and the group 

r = 7ri(E)/7ri(£*) (4.2) 

is the group of covering transformations of £*, interchanging points with the same projections to E. For our 
choice of the covering, T is the free group with TV generators. Choose the generators to be the classes of Bi 
and denote the corresponding transformations of E* by ji. 

This covering arises in the Schottky construction (see Appendix A). In that case T is the Schottky 
group, E* is the Riemann sphere without fixed points of T, 7, are the projective transformations generating 
r. We should emphasize though, that the construction of this covering does not restrict our choice of 
parametrizations of the surface and the moduli nor the choice of projective structures. 

After introducing the covering E* we can define the twisted 1-form as a 1-form on the covering obeying 
proper transformation rules. Assume that the local coordinates on different sheets of the covering are related 
by the action of T. Then we call f a (z) a twisted 1-form in z, a if for alH = 1, . . . , N 

TiWrMz)) = {9i) a b f h {z), (4.3) 



where (gi) b are the matrix elements of the twists in the adjoint representation. 
Main examples of meromorphic twisted 1-forms. 

To reconstruct twisted 1-forms from their singularities and integrals over A-cycles, we shall decompose 
them in a sum of suitable finite-dimensional spaces. For our construction we need to fix an auxiliary point 
wq on the surface. 

(i) Twisted 1-forms wf b (z;wo)- Consider first the space of all meromorphic twisted 1-forms with the 
only simple pole at wo- This space is (N dim C/)-dimensional, and any 1-form from this space is uniquely 
determined by its integrals over A-cycles. Define a basis w-^z; wo) in this space by 

i uj? b (z;wo)dz = 5%5 ij . (4.4) 

The residues at z = wq are 

Res Wo Lot h (z;wo)^(gr 1 -l) a h . (4.5) 
In the Schottky parametrization such twisted 1-forms are given by the Poincare series 

cot b (z;w ) = J2(gy 1 )b 

(ii) Twisted 1-forms 8 b (z;w,w )- Another useful example of twisted 1-forms are 1-forms 6 b (z;w,w ) 
with zero integrals over A-cycles and a given simple pole at a given point w (which is compensated by a 
simple pole at the auxiliary point wq): 

£ 8Z(z;w,w )dz = 0, (4.7) 

Res w 6%(z;w,w ) = -Res Wo 6%(z; w, w ) = 6%. (4.8) 
In the Schottky parametrization the Poincare series for 6 b (z;w,w ) is 

e a b (z-w 1 w )^Y,^ 1 ) a b 

7er 



Y(z) 



■y(z) - ■jAwq) 7(2) -tu 



(4.6) 



V(*) 
7 (z) 



i{z) 



w 



l{z) - w 



(4.9) 
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Here we assumed that w and wo belong to the fundamental domain of T. Holomorphically extending 
these expressions to other sheets of the covering £*, we can write 

wl{z;w )=6 a b {z; li {w),w ). (4.10) 

The transformation properties of b (z; w, w ) as a function of the parameters w or wq are described by 

6Z(z nj (w),w Q ) =6 a c (z;w,w )(gr 1 y b +u;%(z;w ). (4.11) 

The proof of this equality follows from the observation that l.h.s. and r.h.s. are both twisted 1-forms in z, a 
and have equal singularities and equal integrals over A-cycles. 

(Hi) Twisted 1-forms £l ab (z,w). Differentiating 9 b (z;w, w ) with respect to the parameters we obtain 
another twisted 1-form Q ab (z,w) with a pole of second order: 

n ab (z,w) = —9%(z;w,w ) = -—6%(z;wo,w), (4.12) 

The first order residue of il ab (z, w) at z — w vanishes as well as integrals over A-cycles. A remarkable 
property of Q ab (z, w) is its symmetry under interchanging (z, a) <-> (w,b): 

n ab (z,w) = n ba (w, z ), (4.13) 

in particular, fl ab (z, w) is a twisted 1-form both in z, a and in w, b. 
The Poincare series for Q ab (z,w) in Schottky parametrization is 

W b (zM = ^r b { ^ wr (4-14) 
The twisted 1-forms u)f b (z; wo) and 9 b (z; w, Wo) can be expressed in terms of Q ab (z, w) as 

r~ti{wa) 



u>° b (z;w )= / d£n ab (z,Z), (4.15) 

J Wq 

6 a b (z;w,w )= / dtn ab (z,Z). (4.16) 

J Wo 

In our discussion we shall need to look at these forms at their singular points. For this reason we wish 
to treat the twisted 1-forms we defined above as correlation functions in the twisted b-c system of spin 1. 
Then the regularization of singularities will correspond to the normal ordering of the fields. 

Twisted b-c system. 

The b-c system of spin 1 contains anticommuting fields b of spin 1 and c of spin (i.e. b field transforms 
as a 1-form and c field — as a function on the surface). For constructing the twisted version we take dimQ 
copies of b-c systems. The fields b a and c a obey the OPE 

b a {z)c b {w) = J— + (1), (4.17) 
z — w 

and the stress-energy tensor is given by 

T(z) = - :b a {z)dc a (z) : = - lim ( b a (z)dc a (w) - j^^ ] . (4.18) 

w^z \ (z — W) z ) 

We define b-c currents by 

f{z) - r bc b b (z)c c (z), (4.19) 
6 



Then these currents form the affine algebra Q at the level 2/i*: 

sab £abc 

f{z)j\w) = -lh*j—^ j^fH + 0(1). (4.20) 

On the handles of the surface we insert the twists 

g t = exp^tff(z)dz^ (4.21) 

by the same group elements as in the WZNW model. 

To eliminate zero modes of b and c fields we include in correlation functions, along with the twists, the 
product of all c fields in the auxiliary point wq and the integrals of all b fields over all A-cycles: 

dim g N 

n ca ( w o)ir ? ba ( z ) dz - ( 4 - 22 ) 

a=l i=l J A i 

Similarly to the twisted WZNW model, we shall suppress in notation both the twists (4.21) and the 

insertions (4.22), thus writing (Xf- C instead of (X l\f™ S c a (w ) § Ai h<1 ( z ) dz ll^i 9i) h ~ c - 

Correlation functions. 

Consider the normalized correlation function {b a (z)c b '(u>))^~ c . The twists (4.21) make it a twisted 1-form 
in z, a. Its integrals over A-cycles vanish due to the insertion of J\ i a § A b a (z) dz. It has two poles: at z = w 

Res w (b a (z)c b (w)) b N c = 8 ab (4.23) 

and at z — w (due to the insertion of J[ a c a (w )). Therefore, 

(b a {z)c b {w)) b N c - 9Uz;w,w ). (4.24) 

Also 

n ab (z, w) = (b a (z)dc b (w)) b N c = (b b (w)dc a (z)) b N c . (4.25) 
Many-point correlation functions obey the Wick rule. For 4-point correlation functions it reads: 

(6 a (z 1 )& b (^) C c (^3)c £i (z4)) fc N c -(6 a (^i)c d (z4)) fc N c (b b (z 2 )c c (z 3 )) b N c 

) b N c {b b (z 2 )c d (^ , 



(b a ( Zl y(z 3 )) b N - c (b b ( Z2 )c d (z 4 )) b r 



(4.26) 



This identity can be checked easily by comparing the transformation properties and the singularities of the 
l.h.s. and the r.h.s. We shall need this formula in further computations. 
Regular ization. 

We can regularize singularities by normally ordering the fields inside the correlation functions. Mathe- 
matically, this amounts to discarding the negative power terms in Laurent series, i.e. 

Q ab (z,z) reg = (:b a (z)dc b (z):)» N c = lim (n ab (z,w)- — ) , (4.27) 

e a b (z;z,w ) reg = (:b a (z)c b (z):) b N c = lim ( 6 a b (z;w,w ) - -^-) . (4.28) 

This way of regularization depends on the choice of the local coordinate. Although we shall admit such a 
dependence in intermediate calculations, we can check afterwards that the final result transforms properly 
under changes of coordinates. We shall indicate the regularization by the subscript "reg" . 

Stress- energy tensor and partition function. 

Using this notation, 

{T(z))» N c = -n aa (z,z) reg . (4.29) 



If we choose a proper family of projective structures on the surface (e.g. the projective structure of the 
Schottky representation), the stress-energy tensor can be integrated (on a covering of the moduli space) to 
a holomorphic partition function Z b ~ c [6,5]: 

d m (z)logZ»- c = {T(z)) b N c (4.30) 

(see also section 3). The square root II of this partition function will play an important role in our discussion: 

Z b - C = IT 2 . (4.31) 

Notice that II depends both on the moduli of the surface and on the twists </j. 
The following property of II follows immediately from its definition: 

d m (z)Tl = -^im aa (z,z) reg (4.32) 
In the Schottky representation II can be computed explicitly [6] : 

oo 

nSchottky = T[ TT dct(l- 57 i^). (4.33) 

adj 

k—1 *y£prim 

In (4.33) the product is taken over the primitive conjugate classes of the Schottky group (7 and 7 _1 are 
elements of the same primitive class), the determinants are computed in the adjoint representation, K 7 is 
the multiplier of the transformation 7 (sec Appendix A). However, we shall not further need any explicit 
expression for II, and reproduce the formula (4.33) only to simplify comparing with other works. 



5. Knizhnik-Zamolodchikov-Bcrnard equations. 



KZB equation. 

Arising from the Sugawara construction (2.4), the KZB equation must have the form 



(d m + A) F = 0, 



(5.1) 



where d m is the differential on the moduli space, the connection form A is an operator acting on the twisted 
conformal block F. It seems useful to write equations for the product (FU) instead of F, since in this case 
the dependence on the level k reduces to the factor (k + h*)^ 1 in front of the connection form. 

Let F — (<I>i(zi) . . . $„(z„)) be a non-normalized correlation function (conformal block), with the twists 
implicitly included, II be the square root of the partition function of the twisted b-c system (see section 4). 
Let 9% (z; w, wq), u? b (z; wo) and £l ab (z, w) be the twisted 1-forms defined in section 4, d m (z) be the differential 
on the moduli space which maps functions on the moduli to 2-diffcrentials on the surface (see section 3). 
Then the KZB equation looks like follows: 



d m (z) 



1 



2(k + h*) 



n N n N 

(E Z«> E cia " b ia{z\ «*>)) (E fa ™o)t c + E w °) £jc ) 

=1 



13=1 



3 = 1 



+ U(z) 



j ™ = 0, 



(5.2) 

where all right-invariant derivatives with respect to the twists C la must be thought of as acting also outside 
the brackets (in particular, on FII), and the "potential" term 



U(z) = h*fl aa (z,z) reg - ^C ia u; b ia (z;w )u; b jc (z;w )^ c U 
9h* / 1 \ 

= -—(d m {z) + -£»^(^o)^ c (^o)£ j jn 



(5.3) 
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is just a scalar factor. The derivation of this equation is presented in Appendix B. To make the structure of 
the equation more transparent we need to introduce more notation. Namely, we generalize the construction 
of pairing 1-forms with O-boundaries to the twisted setup at higher genera. 
Pairing 1-forms with O-boundaries. 

First recall how one can pair 1-forms with O-boundaries. Let w be a holomorphic 1-form on a simply 
connected domain. Let C be a 0-boundary on this domain with values in a linear space L, i.e. a set of points 
z\, . . . ,z n labelled with elements l\, . . . , l n of L such that 



Define the pairing between C and lu 



h + ... + l n = 0. (5.4) 
uj, (5.5) 



/ "• P2 



taking values in L. Here wq is an arbitrary auxiliary point inside the domain. Condition (5.4) along with 
the holomorphicity of u> makes this definition independent of the choice of wo and pathes of integration. In 
other words, we integrate uu over a 1-chain with the boundary C, and the result does not depend on the 
particular choice of the 1-chain. 

This definition obviously extends to the case of multiply connected domains and meromorphic forms 
u>, if all integrals of oj along noncontractible contours and first order residues at all poles vanish. We shall 
further extend this construction to twisted 1-forms on Riemann surfaces, but first let us look more closely 
at the KZ equation to motivate our discussion. 

Usually we write the KZ equation in the form (1.1) although this expression is not quite invariant under 
the action of SL(2, C) projective transformations. Without using invariant properties of the conformal block 
F = ($i(zi) . . . Q n {z n )) equation (1.1) gets transformed to the most general form 

f3=£a / 

Thus the equation (1.1) uses the choice of the auxiliary point wo = oo. Of course, equations (5.6) and (1.1) 
are equivalent on the class of (/-invariant functions 




a=l 



therefore we may choose any convenient value for wq. The proof of (5.7) is obvious: we surround each point 
z a by a contour so that the sum of these contours is homologically equivalent to zero. Then, integrating the 
current j a (z) along these contours we obtain 

n n . , 

(E*«) F =<£f ^-f{z)$ 1 {z 1 )...3> n {z n ))=Q. (5.8) 

a=l a =l^ Za 

To make the independence of wq explicit in notation we rewrite the equation in terms of pairing with 0- 
boundaries. Namely, the points z a labelled with the operators form a 0-boundary with values in the space 
of operators on (/-invariant elements / £ V\ <g> . . . <g> V n , where V a is the representation space of the field <I> Q . 
We denote such a 0-boundary by C a . Then the KZ equations take the form 



A + r dz 

8z a k + h* a \J {z a -zf 
c° 



reg 



F = 0. (5.9) 



Here dz/(z a — z) 2 is a meromorphic 1-form on the sphere with the only pole of second order at z — z a and 
with zero first order residue. The subscript "reg" means discarding singular powers when we integrate to 
the singular point z a : 

dz \ / f dz ft \ ,„ „„. 

[ .^ --rM- (5-10) 



(z a -z) 2 Jreg z' a ^z a \J {z' a - z) 2 

c c a 
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The rcgularization depends on the choice of the coordinate z around z a , which is consistent with the trans- 
formation properties of <f> a as a conformal field with the conformal dimension 

j-a j-a 

A - = W?F-y (5 - n) 

Generalization to higher genera. 

This notation admits a generalization to higher genera. Now the consideration analogous to (5.7) leads 
to the C/-invariance condition 

n N 

(E £ - E( £W - K a ))F{9) = 0. (5.12) 

QL — 1 i — 1 

Here F is the twisted conformal block depending on the twists gi, N is the genus of the surface. C ia and 
lZ la = (g^~ )%C lb are the right- and the left-invariant derivatives along the twists ((g^ 1 )^ arc the matrix 
elements in the adjoint representation). The analogue of the 0-boundary of the previous construction will be 
a new object ("generalized 0-boundary") consisting of the points z a labelled with t° and the operators C la 
labelling handles. Again, denote this object by C a . If u a is a meromorphic twisted 1-form with vanishing 
first order residues and zero integrals over A-cycles, then the operator 

^ a = E/ w X+E/ uja£ia ' ( 5 - 13 ) 

a=1 Jwo i=l J ™o 

when restricted to (/-invariant functions (5.12) is independent of the choice of wq and paths of integration 
(summation over a is assumed). The integrals in (5.13) are defined on the covering £*, therefore there is no 
ambiguity in the projections of the homology classes of integration pathes onto B-cycles. The arrow above 
the integral sign indicates that we place differentiation operators C la to the right, thus they do not act on u 
itself. Similarly 

^ a = Ew ^ a + E £W / ^ ( 5 - 14 ) 

a=l J w ° i=l J w ° 

is the conjugate (up to sign, since t a a and C la are antisymmetric) operator. 

More formally, in our construction we integrate 1-forms over an element of the operator-valued T- 
invariant relative homology group H{ ' (T,* , M, L), where S* is the "Schottky-type" covering (see section 4), 
M is the set of preimages of the marked points on E*, L is the space of operators acting on twisted conformal 
blocks, r acts on 1-chains by mapping in E* and simultaneously conjugating by the corresponding twist g 1 . 
The element of h[ T) (E*,M, L) as a 1-chain on S* has the boundary consisting of the marked points with 

(V) 

the corresponding operators. This condition determines the cycle in H\ (£*, M, L) unambiguously after we 
fix the basis of A- and B-cycles on the surface by placing twists on the handles and specify the B-projections 
of the cycle to be the differentiation operators along the twists*. 
Final form of the equation. 

Using this notation we rewrite the KZB equations as 

(d m (z) + 2{k l +hst) [ J n ab (v, z) dn J n bc (z, d( + U(z)] J (FU) = 0. (5.15) 

Now recall that Q ab (z,w) = (b a (z)dc h {w)) b N c in the twisted b-c system. Therefore, we naturally define 

d m (z)n ab (u, w) = (T(z)b a (u)dc b (w)) b N c - (T(z)) N (b a (u)dc b (w)) b N c . (5.16) 



* The author thanks A.Losev for elucidating this point. 



10 



Since T = 



— :bdc:, the Wick rule (4.26) gives 



d m (z)n ab (u,w) = n ac (u,z)n cb (z,w). 



(5.17) 



This enables us to rewrite (5.15) as 

1 



d m (z) + 



2(k + h*) 



d m {z)(j d V J d(n ab ( v ,0) re ^j +U(z) 



C C"> 



(FU) = 0. 



(5.18) 



Here the regularization is introduced to eliminate logarithmic divergences when we perform both integrations 
to the same marked point z a : 



dn j d^ ab (r,,0) re =/™[f dv f d^ ab { V ,t) - Y^t a a t a a \og{z' a ~ z a ) 



a=l 



(5.19) 



Recall that a connection V = d + \A is flat (the equations (d^ + XA^)F = are compatible) for any A if 
and only if dA = and [A^, A v \ = for any directions \i and v on the moduli space. The first of these two 
conditions is almost checked (the first summand is locally integrated), except the potential term U(z). U{z) 
is a holomorphic 2-differential on the surface, and therefore represents a 1-form on the moduli space (since 
it can be coupled to Beltrami differentials, see section 3). Appendix C contains the proof that this 1-form 
on the moduli space is closed, i.e. d m (w)U{z) — d m (z)U(w) — 0. Therefore, there locally exists a function 
W of the moduli and the twists such that 



U(z) = d m (z)W. 



(5.20) 



We may define the "universal" operator 



ikzb 



dr] J d£ft 

c b 



ab 



w. 



(5.21) 



Using this operator, we rewrite the KZB equations in the most invariant form 



d m + 



1 



2(k + h*) 



) (Ffl) = 



(5.22) 



for any variation of the moduli. 

Dependence on the coordinates. 

Let us discuss how the terms of this equation depend on the choice of the coordinates on the surface. 
Notice that d m AKZB is not quite a 1-form on the moduli space (or, equivalently, Akzb is not a function). 
Indeed, there are two coordinate dependent regularizations in this term. One of them is performed in (5.19) 
and depends on the local coordinates at the marked points. This dependence gives the correct conformal 
dimensions (5.11) of the fields $ Q . The other regularization appears in U(z) and depends on the projective 
structure on the surface. This regularization provides for d m A K zB the transformation properties (3.4) of 
the stress-energy tensor with the Virasoro central charge 



h* dim Q k dim Q 



k + h* 



k + h* 



dim Q 



"wZNW 2° b 



(5.23) 



which is exactly the central charge of (-FII) . 

Remark that since U(z) depends on the choice of the projective structure, the question of whether it is 
zero does not make sense. However, if U(z) does not depend on the twists, then there exists a projective 
structure (induced by U(z j) in which U(z) = 0. This happens in the case of torus, but we believe that at 
higher genera U(z) depends on the twists nontrivially. 



11 



Examples. 

We can illustrate our construction by three simple examples. In the case of zero genus (sphere) we have 

n ab (z, w) sphere = 6 (5.24) 
(z — wy 

and 

ATzb=Y.W°^~ z p)- ( 5 - 25 ) 
In the case of torus (genus one) the potential term vanishes in the Schottky parametrization[3]: 

U&hZtky = 0. (5.26) 

If there are no marked points, a torus is characterized by a single complex moduli parameter q. The torus 
with a given q is defined as the quotient of C* by the equivalence z <~ qz. Elementary computations show 
that on a torus without marked points 

A% r z u £=logqA, (5.27) 

where A = C a C a is the Laplacian on the group G of the twist [3]. 

In the abelian case (G = U(lj) U docs not depend on the twists, and the potential U(z) vanishes. In 
this case 

A°kzb = £> l CP ^ (5.28) 

where is the period matrix of the surface: 



pi 

J Wo 



U3. 



l (z)dz (5.29) 



(in the abelian case u)°" b (z; wo) becomes proportional to 5 b and independent of w ). 
Symmetricity. 

Our computations are consistent with Losev's observation that the KZB connection form must be 
symmetric [5]. The symmetricity of Akzb and, therefore, of the connection form d m AKZB indicates that 
the connection preserves a certain pairing between solutions of KZB equations with central charges k and 
(-k-2h*). 

Let solutions F k and F^k-2h' of KZB equations for central charges k and (—k — 2h*) respectively take 
values in tensor products V\ ® ... ® V n and V£ (g> . . . <S> V* of representations of G (here V* is the representation 
dual to V a ). Then we can define a pairing between Fk and by 

{F k ,F^ 2h ,) b - c = J dgn 2 (g)(F k ,F_ k _ 2h *) 

, n (5.30) 



(AO 

N3N ' 

1=1 



J] d 9i n 2 (.g)^— (g)F^ h , (g^i --4 



where rj^ are the matrices of the pairings V a ® V* — *■ C, the integration is performed over all twists, dgi is 
the invariant measure on G. With respect to this bilinear form operators t%, , D a and 1Z la are antisymmetric, 
therefore Akzb and d m AKZB are symmetric. Hence, for the solutions F k and F- k -2h' of the KZB equations 
we have 



d m (Fk, F-k-2h' ) b c — (d m F k , f 1 _/ c _2/i*)'' c + (Fk, d rn F- k -2h* ) b c 

1 / \b-c 1 / \b-c (5.31) 

= ]^{( d mAKZB)F k ,F_ k _ 2h ,) + { _ k _ 2h , ) + h , [F k , (d m A KZ B)F- k - 2h ,) = 

due to the symmetricity of d m AKZB- Losev [5] interpreted the symmetricity of the connection as a conse- 
quence of an operator formalism in G/G WZNW coset theory. Referring to the operator formalism we could 
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treat F_k-2h* as a conformal block with the central charge (— k — 2h*), and the twisted b-c system — as 
ghosts in BRST construction [5]. 
Integrability. 

Now let us return to the integrability conditions. To explicitly prove the flatness of the connection it 
remains to show that 

\d il A KZB ,d v A KZB \={S (5.32) 

for any directions \i and v in the moduli space. In Appendix D we prove this for <9 P and d v being derivatives 
with respect to the coordinates of inserted fields. We also derive the corresponding generalization of the 
classical Yang-Baxter equation, which though does not seem to be very instructive (see Appendix D). Instead 
of commuting the most general operators with marked points, in Appendix E we prove the compatibility 
of the equations for the partition function (with no marked points). Wc claim that since marked points 
can be obtained in a degeneration of the surface (double points), we indeed proved the general form of 
the statement. The compatibility of the equations for the setup with marked points (and the result of 
Appendix D in particular) can be obtained as a limiting case of the statement of Appendix E (we present 
an independent derivation in Appendix D just to study the generalization of the classical Yang-Baxter 
equation). Thus we proved that the KZB connection is flat. Remark that our proof did not refer to any 
particular projective structure; the only property we used was that the holomorphic b-c partition function 
exists (i.e. the projective structure is "compatible" — see section 3). 

6. Conclusion. 



To summarize, we have written the KZB equations and proved their integrability. The concise form 
(5.21), (5.22) of the equations contrasts sharply with our way of proving the compatibility (Appendices 
C,D,E). We still hope that the structure of the equations and their relation to the b-c systems can suggest 
a more elegant treatment of the problem and help to avoid the tedious computations. Another question 
remains open about the meaning of the potential term. One may look for its geometric description or for 
its expression in the WZNW coset construction [5] . We also admit the possibility to explicitly integrate the 
potential U(z) to W (see equation (5.20)). Another approach to the problem can be inspired by the works 
on geometric quantization of the moduli space of flat connections, which arrive to similar equations [8]. 
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Appendix A. Schottky parametrization of Ricmann surfaces. 

In the Schottky parametrization the surface is constructed as the quotient of the Riemann sphere (more 
strictly, of the sphere without the fixed points of the group) by the action of a Schottky group. The Schottky 
group r is a group freely generated by N projective maps 7$ such that one can find 27V circles A4 and 
A\ = -fi(Ai), i = 1, . . . , N — all external to each other, and 7, maps the exterior of Ai onto the interior of 
A\. The exterior to all the circles Ai and A\ is a fundamental domain of T. The surface is obtained by gluing 
each circle A\ to Ai by the action of 7$ . Then Ai become A-cycles on the surface. 

For future use we should introduce two more definitions. The first one is the parametrization of a 
projective transformation 7 by its fixed points w 7 , v 1 (repulsive and attractive) and its multiplier K 1 defined 

by 

44^2 =K, Z —^, \K 7 \<1. (AI) 
7(z) — v 7 z — v 7 

Finally, we shall need to extend the twists introduced on the handles of the surface to the group 
homomorphism between T and G: 

S : T - G, 

9ji = 9i, (A-2) 
9ip. = 9i9v 
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The global coordinate z on the Ricmann sphere defines naturally a projective structure on the surface. 
It is remarkable that in this family of projective structures the stress-energy tensor of b-c systems and of 
the WZNW model can be integrated to partition functions (or conformal blocks) according to (3.3); i.e. in 
Schottky parametrization the KZB connection becomes flat. 

Appendix B. Derivation of the KZB equations. 

To derive the KZB equations we first express correlation functions with inserted currents in terms of 
the twisted 1-forms 6 b (z;w,Wo), ujf b (z;wo), n ab (z,w), introduced in section 4. Let $ denote the product 
Then 

a i 

since r.h.s. and l.h.s. are both twisted 1-forms in z, a, have equal singularities at the poles z = £ a and equal 
integrals over A-cycles. Although each of the two terms on the r.h.s. depends on the choice of the auxiliary 
point wq and has a pole at z = wo, their sum does not. For further convenience we introduce more notation: 

L a (z;w a ) =LU* b (z;w )C b , (B.2) 

L a (z; w ) = C ib u? b (z; w ) = L a (z; w ) + {£ ib u,? b (z; w )}. (B.3) 

Here we use braces to specify the range of action for the differentiation operators C la . By convention, all 
differentiation operators act only inside the braces. Also the summation over all repeating indices (including 
those labelling marked points and handles) is assumed if not specified otherwise. 
For the correlation function with the insertion of two currents we have 



<j°(*)j»*> = - kn ab (z,wm e a d (z;w,w )f Mc (j c (w)$) 

+ e a c (z; Z a ,w Q )t c a (j b (w)$) + L a {z- w Q ){j b (w)^) 



for the same reason of the equality of singularities and integrals over A-cycles. Using (B.l) this becomes 

(j a (z)j\w)$) = -kQ ab (z, w){$) 

-e a d (z;w,w )f dbc e c e (w;Z a ,w K(<S>) 
-e a d (z;w,w Q )f dbc L c (w;w )(<S>) 

+ 9 a c (z; £ a ,w a )9 b d (w; (5.5) 
+ a c (z;^w Q )t c a L b (w;w Q )(<S>) 
+ L a (z;w )9 b c (w;^,wo)t c ^) 
+ L a (z;w )L b (w;w Q )(^). 



Introduce 

Then, taking the limit w 



~ a (z;w ) = f abc e b c (z;z,wo) = (j a (z)) b ^ c . (5.6) 

(:j a (z)f(z):<l>) =-kn aa (z,z) reg ($) 

+ E a (z;wo)9 a e (z^ a ,woK(^) 
+ Z a (z;w )L a (w;w Q )(<I>) 

+ 9 a c (z; w )6 a d (z; w )t c a t d ($) {B.l) 
+ e a c (z^ a ,w )t c a L a (z;w a ){<f) 

+ L a (z;w Q )9 a c (z;^,w )tf 3 ($) 
+ L a (z;w a )L a (z;w )(<I>). 
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Now we use Losev's observation [5] that 

~ a (z; w ) = {2L a (z; w ) logH + C b Lu? b (z; w )}. 



(B.8) 



The proof is based on checking that r.h.s. and l.h.s. are twisted 1-forms with additive B-periods with equal 
jumps at A-cycles and integrals over A-cycles. 
This allows us to rewrite (B.7) as 



(:f(z)j a (z):$)= [-kn aa (z,z) reg 

+ (l a (z- w ) + t b a a b (z; fa w )) (l°(z; w ) + t c 6 a c (z; fa w j) 
+ 2{L a (z; w ) logn} (L a (z; w ) + t c p 9 a c (z; fa woj)] <$). 
By the Sugawara construction (2.4) and the identity (4.32) we have 

n(Z Q (z; w ) + t h J a b {z; fa w )) (L a (z; w ) + t c p 6 a c (z; fa w j) 



.(*)(<*>n) 



+ 2{L a (z; w )U}^L a (z; w ) + ^(z; fa w j)] ($> 

-^y([^(^ W 

+ (l a {z- wo) + t h J a b (z; faw )) (L a (z; w ) + t%d a c {z- fa w ))] (<$)n) 
- ($}L a (z;w )L a (z;wo)IlJ 

^— ( (Z a (z; wo) + t b J a b {z- fa w )) (L a {z; w ) + t c 9 a c (z; fa w )) + U(z) j (($)n) , 



2(k + h*) 



where 



U(z) = {h*n aa (z, z) reg - -L a (z; wo)L a (z; w )Il}. 



Appendix C. "Potential" term is a closed 1-form 
on the moduli space. 



Closeness of the potential term U (z) is equivalent to 

d m (w)U(z) - d m (z)U(w) = 0, 



(B.9) 



(5.10) 



(£•11) 



(C.l) 



where the operators d m (z) should be understood as in (5.16), (5.17). Before proving this let us derive several 
useful identities. We shall use the notation introduced in Appendix B. 
First, we wish to prove that 



L a (z; w )n bc (w, - L\w- w )n ac (z, £) 
+ \n ad (z, i)6 b e {w- C wo) - n bd (w, $6 a e (z; £, wo) 



f 



cde 



(C.2) 



n ec (w, £,)0 a d {z; w, w )f bde - tl ec (z, O0 b d (w, z, wo)f aa 







The proof consists of the following steps: (i) Observe that l.h.s. is a twisted 1-form in z, a. To check this we 
use (4.11) along with the identity 

gt^g- 1 = (g) b a t b (C.3) 



15 



in the adjoint representation, (ii) L.h.s. is regular at z — w. (iii) L.h.s. is regular at z = £. (iv) Integrals 
over A-cycles vanish. Assuming that u>, £, wq belong to the fundamental domain, 



(l.h.s.)dz = C la fl bc {w, f ) - (f> fl ec {z, Z)9 b d (w; z, w )f aae dz 



A i 



ade , 



Ai 



(C.4) 



We can prove that this is identically zero by considering correlation functions in the twisted b-c system: 

(C.5) 



c a n bc (w,o = £ ta (b b Hdc c (0) b N 



( d> f(z) dz b»(w)dc c (0) b x c ( f f (*) dz)° N c (b»(w)dc c (0) D x 



\b-c iubi 



\b-c 



Ai 



Recalling that j a (z) = f abc b b (z)c c (z) and using the Wick rule, we arrive to 

C a n bc (w,0 = -f ade j A dz{b b (w)cZ(z)) b N c (b d (z)dc c (t)) b N c 



f ade f A dzn dc (z,oo b e (w^^o), 



(C.6) 



which finishes the proof. 

We can obtain more identities by integrating (C.2) in £ along different pathes. If £ runs from w to w, 
we have 

L a (z; w )9 b c (w; u, w ) - L b (w; w )9^(z; u, w ) 

+ 8 a d (z;u,w )6 b e (w;u, Wo )f de (C.7) 
+ 9 e c (w; u, w )6 a d (z; w, w a )f bde - 9 e c {z- u, w )9 b d (w; z, w )f ade = 



Setting u = ^(wo), this easily leads to the commutation relation 

[L a (z; w ), L b (w; w )} = f acd 9 b c (w; z, w Q )L d (z; w ) - f bcd 9 a c (z; w, w Q )L d (w; w ). 



(C8) 



Let us now return back to proving (C.l). Since we are interested in the antisymmetric with respect to 
interchanging z <-> w part, we shall use the symbol ~ to indicate that the antisymmetric parts of the two 
expressions are equal. 

d m (z)U(w) = h*d m (z)Q aa (w, w) reg - ^d m (z) (^L a (w; w )L a (w; w )I^j . (C.9) 

The first summand has zero antisymmetric part, since Q aa (w, w) reg — —2d m (w) log II is an exact 1-form on 
the moduli space, and d^ = 0. Also, since 



, pi(w ) pj(wa) . 

L a (w;w a )L a (w; Wo ) = d m ( W )(C b dU dry O bc (£, t?)/}^) 

J Wn J W(\ 



(CIO) 



we arrive to 
d m (z)U(w) 



(d m {z)^)L b (w : w )L b {w; w )n + ^L b {w; w )L b {w; w ) (d m (z)Il) 



_ ^ y_LJl51lb {w . WQ ) L \ W . WQ )U + -L b (w; w )L b (w; w ) (im aa (z, z) r 

i [^{i V w Q )n aa (z, z) reg }{L b (w; w )IL} + L b (w; w Q )L b (w; w )n aa (z, z) 
X - \L b (w; w )L b (w; w )n aa (z, z) reg + E b (w; w )L b {w; w )n aa (z, z) 



(C.ll) 



)reg 
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From (C.2) we deduce that 

L b (w; w a )n aa (z, z) reg = L a (z; w )n ba (w, z) + [n ea (w, z)9 a d (z; w, w )f bde + fl bd (w, z)E d (z; w ) 
Using it in rewriting the first summand in (C.ll), we obtain 

4d m (z)U(w) ~ L b (w;w )L a (z;w )n ab (z,w) 

+ L b (w; w ) [n ea (w, z)9 a d (z; w, w )f bde + n bd (w, z)E d (z; w )] ~ E b (z; w )L b (z; w )n aa (w, w) reg 

Commuting L a (z;w ) and L b (w : wo) in the first term according to (C.8), 

Ad m (z)U{w) ~ f bcd 9 a c (z; w, w Q )L d (w; w Q )n ab (z, w) - f edb 9 a d (z; w, w )L b (w; w )(l ae (z, w) 
+ n ea (w, z)L b (w; w )9 a d (z; w, w a )f bde 
+ n bd (w,z)L b (w;wo)E d (z;w ) 



(C.12) 



(CIS) 



E b (z; w ) [L b (z; w a )n aa (w, w) reg - L a (w; w )n ba (z, w) 
= n ea (w, z)L»(w; w a )8 a d (z; w, w a )f bde 
+ n bd (w,z)L b (w;wo)E d (z;w ) 
-E b (z;wo)n ea (z,w)9 a d (w;z,w a )f bde 
-~ b (z;wo)n bd (z,w)~ d (w;w ). 

The last term is symmetric, therefore, 

4d m (z)U(w) ~ n ea (w,z)L b (w;w )9 a d (z;w,wo)f bde 
- Q ea (w, z)L a (z; w )6 b (w; w, w )f bde 
-Z b (z;w a )n ea (z,w)8 a d (w;z,wo)f bde 

dfg 



(C.14) 



= n ea (w, z) 9 a f (z; w, w )O b (w; w, w )re g f 
+ 9 9 d (w; w, w )re g O a f (z; w, w Q )f bf9 - e g d (z; w, w )9 b f (w; z, wo)f afg 
-Z b (z;w )n ea (z, W )9 a d (w;z,w )f bde 
~ n ea (w, z)9 a f (z; w, w )9 b (w; w, w )re g f de f df9 

+ n ea (w, Z)9 g d (w- W, Wo)re g e a S {z- W, W ) [-f b f d f b ^ - f*°t ' fbdg 

+ Q ea (w, z)E b (w; wo)9 a d (z- w , w a )f bde 
= 0. 



(C.15) 



We used (C.7) and the Jacobi identity 

jabc jade _|_ jabd jaec _|_ jabe jacd 



o. 



(C.16) 



Appendix D. Compatibility of the KZB equations 
for moving marked points. 

We check the compatibility of the equations for the derivatives of the conformal block with respect to 
the coordinates of marked points. It is slightly simplifies computations if we deal with the equations for 
"bare" conformal blocks F — ($i(£i) • • • 3>n(£n))i not for (TIF). The equations for F look like follows: 



d 

d^ a k + h 



(5 (& (reg) */3 + L a (^ a ] Wofj 



0. 



(DA) 
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The regularization of 6% (£ Q ; £3, wq) is assumed when (3 = a. Our goal is to show that the n-point operators 
commute: 

The commutativity property is sufficient to check only for 2- and 3-point operators. The structure of the 
operators ensures that any n-point operators will then also commute. One easily observes this property 
for the KZ equations (1.1). In this case the commutativity of the 2-point operators is trivial, since 

A {2)s P here = _ A (2)sphere (£) ^ 

for two marked points a and (3. The 3-point operators commute due to the classical Yang-Baxter equation 

lKt ere > R 7i ere } + [ R t ere > R t ere \ + l R i<* here . = 0, (i>.5) 

where the R-matrix is 

< m = T^V- (D-6) 
This argument slightly changes in the case of higher genera. Let 

-^i Q (ea^0)+^(ea;e Q ,W0W^, (I>.7) 

fla/J = (I>.8) 

(no summation over a and (3). Then 

4 n) = 5 a + E ^ (£>-9) 

The commutativity of the 2-point operators implies 

[S a +R a/3 ,Sp + R 0a ]=O. (D.10) 
For the 3-point operators the commutativity means 

[A& + R ai , Af + R 01 ] = 0, (£>.ll) 

(2) (2) 

where A a and A^ are the 2-point operators for the marked points £ Q and Up to (D.10) this is equivalent 
to 

[4 2) , RfJj} + K 7 , A^ 2) ] + [R a7 , R 07 ] = (D.12) 
or being rewritten in another way: 

[R a /3,R/3j] + [R ai ,R 01 \ + [R ai ,R 0a \ = [R/3 7 ,S a ] + [Sp,R a7 ]. (DAS) 

This is a generalization of the classical Yang-Baxter equation (D.5) to higher genera. Notice that here we 
do not have the symmetry R a0 = —Rp a any more. 

Our proof of the compatibility condition (D.3) will proceed in two steps. First we check (D.13), then 
(D.10). The generalized classical Yang-Baxter equation (D.13) follows from (C.7): 

[Ra@, Rp~(] + [Raf, Rp-i] + [Ran Rp a ] ~ [Rp-y, S a ] ~ [Sfj, R al ] 

+ ^(e a ;^, Wo )^(^;e7^o)/ bde -^(^;e Q , W o)^(C a ;e 7 ,wo)r de (£.14) 

+ L a (^;w o )0 b c (^; ^,w ) - L b (faw )6 a c (£ a ;^, w )} 
= 0. 
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Commuting the 2-point operators, obtain 



[S a + R al 3, Sp + Rpa\ = Kf c t^tpC c + X a c tat a tp + X a c t a pt fi t c al (D.15) 



where 



(£.16) 



+ L a {^ a ;w Q )uj b ic (^; w Q ) - L b (&; tuo)wg.(^ Q ; w ) 
-0, 

X abc {^p) = -X ahc (£p,Z a ) 

= f bde 0d (&» 5 £a , W ) reg 6 C e (& ;C„,W0) + / ad6 ^ J £a , l«o)re fl ^ (& ;&,«*>) 

- L c (& ; w ) 0£ (&, ; £ Q , w ) reg + L a (&, ; w )0 g ; £a , w ) (D. 17) 

= 

again due to (C.7). This completes our proof of (D.3). 

Appendix E. Compatibility of the equations 
for the partition function. 

We explicitly show that in the case of no insertions (equations for the partition function) the KZB 
equations are compatible (the connection is flat). This follows from 

[d m (z)A K zB,d m (w)A KZ B] = 0. (£7.1) 

We prove it in this Appendix. 
Let 

A(z) = d m (z)A KZB = A B (z) + U{z), (£7.2) 

where 

A B (z) = L a (z ] w )L a {z;w ) (£7.3) 
is a symmetric second order differential operator, U(z) is the potential term, 

u(z) = h*n aa (z, z) reg - {^A B ( z )n} 

= h*n aa (z, z) reg - {A B (z) loglT} - {L a (z) logn} 2 . 

From now on we shall omit the dependence on the auxiliary point wq in all expressions, so that L a (z) = 
L a (z;w ), etc. Commuting A(z) and A(w), we obtain 

[A{z), A(w)] =[A B (z),A B (w)} + [A B (z), U(w)} - [A B (w), U(z)} 

= {A B (z), A B (w)] + {L a (z)U(w)}L a (z) + L%z){l%z)U{w)) (jB . 5) 

- ^L a (w)U(z)^L a (w) - L a (w)^L a (w)U(z)y 
We wish to prove this to be zero. 
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First compute 

[A B (z), A B (w)} = L a (z)L\w)[L a (z), L h (w)} 

+ {L a (z),L b (w)}L»(w)L a (z) 

+ L a (z)[L a (z),L b (w)]L b (w) 

+ L\w)[L a (z),L b {w)]L a {z) 
=L a (z)L b (w)f acd 6 b c (w; z)L d {z) 
-L a {z)L b {w)f bcd e a c {z- w)L d {w) 
+L d (z)f acd e b c (w; z)L b (w)L a (z) 

-L d (w)f bcd 6 a c {z- w)L b (w)L a {z) ( E - Q ) 
+L a (z)f acd 6 b c (w; z)L d {z)L b (w) 
-L a {z)L d {w)f bcd e a c (z; w)L b {w) 

+L a (z)^L a (z)L b (w) + L d (w)f bcd 6 a c (z; w)^L b {w) 

+L b (w)L d (z)f acd 9 b (w; z)L a {z) 

-L b {w)f bcd 6 a c {z- w)L d (w)L a (z) 

-L b {w)^L a {z)L b {w) + L d (w)f bcd 9 a c {z- w)^L a (z). 

In the above transformations we used the commutation relation (C.8) and the one transposed to it. In the 
last expression the second term cancels the sixth one, the fourth term cancels the ninth one. We denote the 
expression in braces by P ab (z, w). Notice that 

P ab (z, w) = {L a (z)L b (w) + L d (w)f bcd 6 a c (z; w)} 

= {L b (w)L a (z) + L d (z)f acd e b (w; z)} = P ba (w, z). 

Before going further let us state a simple consequence of the identity (C.7). Namely, setting u — > w we 
obtain 

L a (z)e b c (w;w) reg - L b (w)9 a c (z;w) 
= -9 a d (z;w)9 b e (w;w) reg r de 

-e a d (z;w)0 e c (w;w) reg f bde (E.8) 
+ 6 b (w;z)e e c (z; W )r de 

+ f bcd n ad (z,w). 

This lemma allows us to rewrite P ab (z,w) as 

P ab (z,w) ={-2L b (w)L a (z)\ogU- L b (w)f acd 9 d c (z;z) 

+ r d b (w, z)L d (z) + L d (z)r d e b (w; z)} 

={-2L b (w)L a (z)\ogn-2f acd e b c (w;z)L d (z)\ogIl 

- f acd f bfe 9 d f (z;w)8 e c (w;z)+2h*n ab (z,w)}. 

We shall further use this expression, and now return to (E.6), which after cancellations and regrouping 
the remaining terms becomes 

[A B (z), A B (w)} ^[L a (z) 1 L b (w)\r d e b c (w- z)L d (z) - L d (z)f acd 9 b (w; z)[L a (z) , L b («,)] 

+ L a (z)P ab (z, w)L b (w) - L b {w)P ab {z, w)L a (z) 

=Lf(z)rfe b e ( W] z)r cd e b c (w; z)L d (z) 

- L d (z)f aef e b (w; z)f acd 6 b (w; z)L f (z) ( E - 10 ) 

+ L d (z) (P d f{z, w) + f bef 6 a e (z; w)r d 6 b (w; z)) \w) 
+ U{w)(p df {z, w) + f bef 6 a e {z; w)f acd 6 b c (w; zfj L d (z). 
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The first two terms of the last expression cancel each other due to the Jacobi identity (C.16). Introduce now 
one more quantity Q ab (z, w) defined by 

Q ab (z, w) = P ab (z, w) + f adc f bfe 9 d (z; w)9{ (w; z). (£.11) 

Like for P ab (z,w) we have 

Q ab (z,w) = Q ba (w,z). (E.12) 

The expression (E.9) for P ab (z, w) gives 

Q ab (z, w) = 2{h*fl ab (z, w) - L b (w)L a (z) log II - f acd 9 b (w; z)L d {z) logll}. (£.13) 
With this notation the commutator (E.10) becomes 

[A B (z),A B (w)} =L a (z)Q ab (z,w)L b (w) -L b (w)Q ab (z,w)L a (z) 

={L a (z)Q ab (z, w)}l\w) - [l b {w)Q ab {z, w)}L a (z) + Q ab (z, w)[L a (z), L b (w)} 
={L a (z)Q ad (z, w) - Q ab (z, w)f bcd 9 a c (z; w)}l a (w) 
- {L a (w)Q ad (w, z) - Q ab (w, z)f bcd 9 a c (w; z))L d {z). 



(£.14) 



This is a first order differential operator. It is antisymmetric although this is not seen at once from (E.14). 
Therefore we can restrict our attention to its symbol only. Returning to the original problem (E.l), we arrive 
to 

[A{z),A{w)\ =Y d (z, w)L d {z) + L d (z)Y d (z, w) 
-Y d (w,z)L d (w)-L d (w)Y d (w,z), 

where 

Y d (z,w) = ^L d (z)U(w) - ^L a (w)Q ad (w,z) - Q ab (w,z)f bcd 9 a c (w;z)^Y (£.16) 
We shall prove that 

Y d (z,w)=0, (£.17) 

then (E.l) will follow. 

Writing out Y d (z 1 w) explicitly, 

Y d (z,w) =h*{L d (z)n aa (w,w) reg ^L a (w)n da (z,w) + n ab ( W ,z)f bcd 9^(w;z)^ 
-L d (z)^L b (w)L b (w)u) 

+L a (w)L d {z)L a (w)lo g n (£.18) 
+L a (w)f acb 9 d (z; w)L b (w) log II 
-f bcd 9 a c (w;z)L b (z)L a (w)\ogU 
-f bcd 9 a c (w; z)r e f9 b e (z; w)L^ {w) log II. 

Transforming the expression in braces, we arrive to 

{L d (z)fl aa (w, w) reg - L a (w)n da (z, w) + Q ab (w, z)f bcd 9 a c (w; z)} 
= \-n da (z,w)L a (w) + tt ba (z,w)f bcd 9 a c (w;z) 

(£.19) 

+ n ea (z, w)9 a f {w- z)f dfe + fl df (z, w)E f (w)j 
= ^2n df (z,w)L f (w)\ogUy 
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Two next summands in (E.18) become 

L a {w)L d {z)L a {w) logIl-L d (z)^L b (w)L b (w)U^ 

= [L a (w),L d (z)]L a (w) log n-2{i a H logllj L d {z)L a (w) log n 

=L b (w)f acb 9 d (z;w)L a (w)\ogU 
-f dcb 9 a c {w- z)L b (z)L a (w) logn 

-{L d (z)L a (w) + L b (w)f acb ef(z;w)}L a (w)\o g n ( £ - 20 ) 
-{2L d (z)L a (w) lo g n}i a H logn 

= L b (w)f acb e d (z;w)L a (w)logU 

- f dcb 9 a c {w- z)L b (z)L a {w) logn 

- {i d (z)S a H +L b (w)f acb 6 d (z;w)}L a (w) logn. 

The first two terms cancel the fourth and the fifth summands in (E.18) thus leading to 
Y d (z, w) =2h*{tt df (z, w)L f (w) log n} 

-{r fa f cbd e{{w-z)e b e {z-w) ( £. 21) 

+ L d (z)E a (w) + L b (w)f acb e d (z; w)}L a (w) logn. 

Notice that 

L d {z)E a (w) + L b {w)f acb 9 d c {z; w) 

= r b ( y L b (w)0 d (z;w) - L d (z)e b (w;w) reg ) (jB . 22) 
+ 6 d c {z; w)f cba (z b (w) - 2L b (w) log n) . 

The very last term {29 d (z; w)f cba L b (w) log n} vanishes after multiplying by L a (u>)logn. The first term in 
(E.22) can be rewritten according to (E.8), which leads to 

Y d {z,w) = 2h*n df (z,w)L f (w)\ogU 

-{f efa f chd e f c {w-z)e b e {z-w) 

+ 9 d c (z;w)f cba E b (w) 

+ r b [e d f (z;w)9 b e (w;w) reg f^ 

+ 9 d (z;w)8 e c (w;w) reg f bfe 
-6 b (w;z)6 e c (z;w)f d f e 

- f bcf n df (z,w)^L a (w) logn 

= 0. 

Here the first line cancels the last one; the second line — the sixth one; the third to fifth lines cancel due to 
the Jacobi identity. 
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